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1 Introduction 

In this work we complete our study of the sub-Riemannian problem on the Lie group SH(2) which is the group of 
motions of pseudo Euclidean plane. The work was initiated in [T] where we defined the sub-Riemannian problem. 
The control system comprises two 3-dimensional left invariant vector fields and a 2-dimensional linear control 
vector. We applied PMP to the control system and obtained the corresponding Hamiltonian system. In [5] we 
proved the Liouville integrability of the Hamiltonian system. We calculated the Hamiltonian flow such that the 
extremal trajectories were parametrized in terms of Jacobi elliptic functions [T]. Since PMP states only the 
first order optimality conditions, the trajectory resulting from PMP are only potentially optimal called extremal 
trajectories or geodesics. Further analysis based on second order optimality conditions is then needed to segregate 
the optimal trajectories or the minimizing geodesics. It is well known that the candidate optimal trajectories lose 
optimality either at the Maxwell points or at the conjugate points |3|,[1],[5]. Based on the optimality analysis 
one is able to state the time of loss of global optimality known as the cut time. Rigorous techniques for this 
optimality analysis have evolved over the years from research on related sub-Riemannian problems on various 
Lie groups, see e.g., El, 0, 0, 0- These techniques were employed in and [5] to compute the Maxwell strata 
and the conjugate locus in the problem under investigation. An effective upper bound on the cut time was also 
computed. 
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In this paper we extend the global optimality analysis similar to [9]. We decompose the image M = SH(2) 
and the preimage of the exponential mapping into open dense sets based on the Maxwell strata and conjugate loci 
and prove that the exponential mapping between these sets is a diffeomorphism. This leads naturally to the proof 
that the cut time is equal to the hrst Maxwell time. Finally, we analyze the global structure of the exponential 
mapping and obtain explicit characterization of the cut locus and the optimal synthesis on the manifold SH(2). 

The paper is organized as follows. In Section 2, we review the results from [1] and [ 8 ] as ready reference. Sec¬ 
tions 3 and 4 contain the main results of this work. In Section 3 we state and prove the conditions for exponential 
mapping being a diffeomorphism and compute the cut time. Section 4 pertains to explicit characterization of the 
Maxwell strata and the cut locus in terms of a stratification of SH(2). In Section 5 we conclude this work. 


2 Previous Work 


2.1 Problem Statement 


The Lie group SH(2) is a 3-dimensional group of roto-translations of the pseudo Euclidean plane |10| . The 
sub-Riemannian problem on the Lie group SH(2) reads as follows [1]; 


i = uicosh2;, y = uisinh2, z = U2, 

? = (r, y, 2 ) G M = SH(2) = a:,y, 2 GlR, (ui,U 2 )gIK^, 
q{0) = (0,0,0), q{ti) = qi = (a;i,yi, 2 i), 



( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


By Cauchy-Schwarz inequality, the sub-Riemannian length functional I minimization problem (2.41 is equivalent 
to the problem of minimizing the following action functional with fixed ti m- 


j = 


(ui -h U2)dt 


(2.5) 


2.2 Known Results 


We now briefly review the results from [I] and [8] as a ready reference in this paper. System satisfies the 
bracket generating condition and is hence globally controllable mm- Existence of optimal trajectories for the 


optimal control problem (2.1|-(2.5l follows from Filippov’s theorem [3]. We applied PMP to (2.1|-(2.5l to 
derive the normal Hamiltonian system. It turns out that the vertical part of the normal Hamiltonian system is a 
double covering of a mathematical pendulum. The normal Hamiltonian system is given as: 


7 = c, c = — sin 7 , A = ( 7 , c) G C = (25.^) x Rc, 2Sj = R/( 47 rZ), 

'V /-y W 

X = cos ^ cosh 2 , y = cos ^ sinh 2 , 2 = sin A. 


( 2 . 6 ) 

(2.7) 


The total energy integral of the pendulum (2.61 is given as: 


E= — -cos^, E G [-1,- 1 - 00 ). 


( 2 . 8 ) 


The initial cylinder of the vertical subsystem is decomposed into the following subsets based upon the pendulum 
energy that correspond to various pendulum trajectories: 


c=[jc,, 


Cl = {A G Cl Eg (-1,1)}, (2.9) 

C 2 = {A G C|E G ( 1 , 00 )}, (2.10) 

C 3 = {AG C|E = l,cyf 0}, (2.11) 

C 4 = {AG C|E = -1, c=0} = {( 7 ,c) G C |7 = 27rn, c = 0}, n G N, (2.12) 

Cs = {A G C|E = 1, c = 0} = {( 7 ,c) G C |7 = 27rn + 7 r, c = 0}, n G N. (2.13) 


where. 













Cut Locus and Optimal Synthesis in Sub-Riemannian Problem on the Lie Group SH(2) 


3 


C 



Fig. 1 Stratification of the Phase Cylinder C of the Pendulum 


We defined elliptic coordinates {(fjk) for A G U^—iCi C C and proved that the flow of the pendulum is 
rectified in these coordinates. Note that k was defined as the reparametrized energy and ip was defined as the 
reparametrized time of motion of the pendulum [T]. Integration of the horizontal subsystem in elliptic coordinates 
follows from integration of the vertical subsystem and the resulting extremal trajectories are parametrized by 
the Jacobi elliptic functions sn{p,k), cn{p,k), dn((p,A:), E{p,k) = dn^{t, k)dt (Theorems 5.1-5.5 ID)- The 
results of integration for A G Ci, i = 1,... ,5, are summarized as: 


— Case I : X = {p, k) G Cl 


/xt\ 

( ^ 

2 

(w+[E((pt) Hv>)] + ku?j[snpt sn:p] 

1 ) 

rr 

1 

2 

(w' »(i-fe 2 )) PCv’t) E(v5)] (^^^y^_^^^+kwj[snpt snvp] 


\ztj 


Si In [(dn:pt — kcnpt).w] 

/ 


where w = dny-\cny ■ = sgn (cos and pt = p + t. 

— Case 2 : A = {ip, k) G C 2 


1 

H 


1 


2 \w(l — 


+ - ( fern + 


-wj |^E(7/)t) - E{ip) - k'^ {ipt - V’)j 
[sn Ipt - sn , 


yt = - - 


S2 


10(1 — k '^) ^ 

sn Ipt - sn V"], 


+ ^ I k'w — 


k 


w(l — 

zt = S 2 In [(dn V>t - kcnipt) .w] , 

where ip=^, ipt = ^ = ip + I &nd w = dnV--fccnV. ’ = sgnc, k' = ^/l-k'C 
— Case 3 : X = {p, k) G C 3 


/ Xt\ 


^ {pt - p) + w (tanh pt - tanh p) 

h/t = 


^{pt- p) -w (tanh pt - tanh p) 


V 

— S 1 S 2 ln[u) sech:/Pi] 


where w = cosh p. 

— Case 4 : A = {p, k) G C 4 


— Case 5 : X = {p,k) G C 5 


'x\ /sgn (cos 5 ) 




^sgn (sin t ^ 


(2.15) 


(2.16) 


(2.17) 


(2.18) 
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The phase portrait of the pendulum admits a discrete group of symmetries G = {Id, ..., The symme¬ 
tries £* are reflections and translations about the coordinates axes ( 7 , c). The reflection symmetries in the phase 
portrait of a standard pendulum are given as: 

: (7,c) -t { 7 ,-c), 

: (7,c) -t (-7,c), 

: (7,c) -t (- 7 ,-c), 

: ( 7 ,c)-t {7 + 27r,c), (2.19) 

-t (7 + 27r,-c), 
e® : (7,c) -t (-7-|-27r,c), 

: (7,c) -t {-7-|-27r,-c). 

According to Proposition 6.3 [T], the action of reflections on endpoints of extremal trajectories can be defined as 
£* : g !->■ 5 *, where q = {x, y, z) G M, q® = (x*, i/*, 2 ®) £ M and, 

{x^,y^,z^) = ( 2 :cosh 2 — ysinh 2 , a;sinh 2 — qcosh 2 , 2 ), 

{x^,y‘^, 2 ^) = (® cosh 2 — 2 /sinh 2 , —xsinh 2 -|- 1 /cosh 2 , — 2 ), 

{x^,y^,z^) = {x, -y, - 2 ), 

{x‘^,y‘^,z^) = {-X, y, -z), (2.20) 

{x^,y^, 2 ®) = (—a; cosh 2 -|- i/sinh 2 , xsinh 2 — 1 /cosh 2 , — 2 ), 

(a;®, 3 /®, 2 ®) = (—a;cosh 2 -|- 3 /sinh 2 , —a:sinh 2 -|- 2 /cosh 2 , 2 ), 

,z^) = {.-X, -y, z). 

These symmetries are exploited to state the general conditions on Maxwell strata in terms of the functions zt 
and Ri{q) given as: 

z z z z 

= 3 /cosh - — a: sinh -, i ?2 = a; cosh - — 3 /sinh -. ( 2 . 21 ) 

We define the Maxwell sets MAX^, i = 1,..., 7, resulting from the reflections £® of the extremals in the 
preimage of the exponential mapping N as: 

MAX® = = {X,t)GN = C X M"*" I A A®, Exp(A,t) = Exp(A®,t)|, 

where A = e®(A). The corresponding Maxwell strata in the image of the exponential mapping are defined as: 


Max® = Exp(MAX®) c M. 

In [H| Proposition 3.7 we proved that the first Maxwell points corresponding to the reflection symmetries of the 
vertical subsystem lie on the plane 2 = 0 and the corresponding Maxwell time t5®^®®®'(A) is given as : 


A G Cl 
A G C2 
A G C 3 U C 4 U C 5 


= 4K(k), 
= AkK{k), 


Similarly we proved that the first conjugate time t““‘'(A) is bounded as (Theorems 4.1-4.3) [ 8 ]: 


( 2 . 22 ) 

(2.23) 

(2.24) 


ag Cl 
Ag C 2 
Ag C4 
A G C 3 U C 5 


4A'(fc) < t“"j(A) < 2p\{k), 
AkK{k) < t‘(°'®j(A) < 2kp\{k), 

= 27V, 

t“"j(A) = +oo. 


(2.25) 

(2.26) 

(2.27) 

(2.28) 


where Pi(fc) is the first positive root of the function /i(p) = cnpE(p) — snpdnp, which is bounded as Pi(fc) G 
{2K{k),3K(k)). Note that we defined: 


pt = r +p, p = T — p 

I ‘fit , I V 

Tpt = — = T + P, -Ip = - = T - P 


X = + , P = ^ when 1 / = (A, t) G NiU N 3 , (2.29) 

r = + = ^ u = (A, t) gN 2 . (2.30) 


Here and below Ni = Ci x R-).. 
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3 Upper Bound on Cut Time 

In this section we describe the basic properties of the upper bound on cut time obtained in [S]. 
Define the following function t : C —^ (0, +oo], 


t(A) = min i A € C. 


Equalities (2.221-( 2.281 yield the explicit representation of this function: 


ag Cl = 

t(A) = 4:K{k), 

(3.1) 

Ag Ca = 

^ t(A) = 4kK{k), 

(3.2) 

Ag C4 = 

^ t(A) = 27r, 

(3.3) 

AGC 3 UC 5 = 

t(A) = -boo. 

(3.4) 


In |H] we proved the upper bound: 

tcut(A) <t(A), AeC. (3.5) 


We now prove that inequality (3.51 is in fact an equality (see Theorem 4.21. The general scheme of the proof 
is as follows m, 0 : 


1. The exponential mapping Exp ■. N = C x R+ —>■ M parametrizes all optimal geodesics, but also all non- 
optimal ones, since all the geodesics Exp(A,t) with t > t(A) are not optimal. 

2. We reduce the domain of the exponential mapping so that it does not include these a priori non-optimal 
geodesics: 

N = {{X,t)€N I t<t(A)}. 

We also reduce the range of the exponential mapping so that it does not contain the initial point for which 
the optimal geodesic is trivial: 

M = M\{qo}. 

Then Exp : N ^ M is surjective, but not injective, due to Maxwell points. 

3. We exclude Maxwell points in the image of Exp: 

M={gGM I e\q)j^q], 

and reduce respectively the preimage of Exp: 

N = Exp"^ . 

The mapping Exp : N ^ M is injective. Moreover, it is non-degenerate since t“"^(A) > t(A). 

4. We take connected components in preimage and image of Exp : 


N = UDi, M = UMi. 


Each of the mappings Exp : Di —> Mi is non-degenerate and proper. Moreover, all Di and Mi are smooth 
3-dimensional manifolds, connected and simply connected. By HjMamard’s global diffeomorphism theorem 
|14| . each Exp : Di —>■ Mi is a diffeomorphism. Thus Exp : —>■ M is a diffeomorphism as well. 

5. Eurther, we consider the action of the exponential mapping on the boundary of the 3-dimensional diffeomor- 
phic domains: 


Exp : N' M', N' = N\N, M' = M\M. 

We construct a stratification in the preimage and the image of Exp : 

N' = \jNi, M' = UM', 
dim Ni dim M' G {0,1,2}, 

where all are disjoint, while some coincide with others. Eurther, we prove that all Exp : —>■ are 

diffeomorphisms by the same argument. 
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6 . On the basis of the global diffeomorphic structure of the exponential mapping thus described, we get the 
following results: 


tcut(A) — A ^ (17, 

Max = U {M/ I 3j^i such that Mj = , 

Cut = cl(Max)\ {qo} , 

Cut n Conj = cl(Max)\ {go} • 

We show that the optimal synthesis is double valued on the Maxwell set Max, and is one valued on M\Max. 
The central notion of our approach is the stratification in the preimage and in the image of Exp : 

N = (UA) U (UM') , 

M = (UM,) U (UM') , 
dim(Di) = dim(Mj) = 3, 
dim(iV'), dim{M') G (0,1, 2} , 

such that all the corresponding strata are diffeomorphic via the exponential mapping, i.e., Exp : Di —>■ Mi 
and Exp : N'i —>■ M[ are diffeomorphisms. 

It is well known mm that for any smooth manifolds X and Y of equal dimensions, a smooth mapping 
/ : Af —> y is a diffeomorphism if /, X and Y satisfy the following conditions PI - P4: 

PI - X is connected, 

P2 - y is connected and simply connected, 

P3 - / is non-degenerate, 

P4 - / is proper, i.e., for any compact set K CY the inverse image f~^{K) C X is also compact. 

We now consider the invariance properties of the function t with respect to the reflections e* G G and the 
vertical part of the Hamiltonian vector field: 

- sin 7 ^ e Vec(C). 
a 7 oc 


Proposition 3.1 

(1) The function t is invariant w.r.t. the reflections e* G G and the flow oflli/: 

t o £*(A) = t o (A) = t(A), A G G, e* G G, f G R. 


(2) Thefunctiont : G —>■ (0,-|-oo] is in fact a function t{E) of the energy E = ^—cosy of pendulum (2.6l. 


Proof The reflections e* G G (2.191 and the flow of preserve the subsets Ci of the cylinder G and on each of 
these subsets, the function t is expressed as a function of the energy E of the pendulum since we have equalities 
(3T|-([0) and. 


A G Gi 

A G G2 

A G G4 
A G Gs U Gs 


E + l 
2 ’ 


k = 

^ ~ '‘E + l’ 
E = -l, 

E = l. 


This proves item (2) of this proposition. Item (1) follows since the energy E Is invariant w.r.t. e* and □ 

A plot of t(E) Is shown in Figure]^ Regularity properties of the function t{E) visible in its plot are proved 
in the following statement. 


Proposition 3.2 


(1) The function t(A) is smooth on Gi U G 2 . 

(2) limB_>_i t(E) = 2n, lim£;_>i t{E) = -|-oo, lim£;_>+oo t{E) = 0. 

(3) The function t : G —>■ (0, -|-oo] is continuous. 
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Proof Item (1) follows from ( |3.1[ ) and ( |3.2[ ). The limits in item (2) follow from ( |3.1[ ) and ( |3.2[ ), and from the limits 
lim/;_>,^Q limfc_^i_ K(k) = +oo. Then continuity of t(A) follows on Ca,: 

A -t A G C4 => E{\) E{\) = -1 t(A) ^ 27r = t(A). 

Continuity on C 3 U C 5 follows since 

A -t A G C 3 U C 5 E{\) E{X) = 1 t(A) -t +00 = t(A). 

Thus t(A) is continuous on C and item (3) is proved. □ 


3.1 Decompositions in the Image of the Exponential Mapping 

Consider the set M = M\{qo}. From Filippov’s theorem and Pontryagin’s Maximum Principle [3], we already 
know that any point q £ M can be joined with go by an optimal trajectory q{s) = Exp(A, s) such that q{t) = 
q, (A,t) G N. Then Exp(A^) 3 M. However the Maxwell points q £ M have non unique preimage under the 
exponential mapping. Hence the mapping Exp : N ^ M is surjective, but not injective. In order to separate 
Maxwell points we consider the set that contains all such points: 

M' = |g G M I z = 0, a;^ + yf o| , 

and its complement M in M: 

M = {q£ M I 2 yf 0}, 

M = MUM', 

where U is the union of disjoint sets. 

3.1.1 Decompositions in M 

The plane 2 = 0 cuts the domain M into two half spaces as: 

M = Ml U M 2 , 

Ml = {q£ M I 2 > 0} , 

M 2 = {qe M I 2 < 0} . 

Note that the decomposition of the manifold M is simpler in description of cut time on SH(2) than similar 
decomposition of M in related problems on SE(2) and on the Engel group [7]. 


(3.6) 

(3.7) 
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Id, 


Ml 

M 2 

M 2 

Ml 


Table 1 Action of e® on Mj 



C 

A 

A 

A 



ti 

1 

0 71 

A 

ItT 

A 

1 


Fig. 3 Projections of Di to Phase Cylinder C of the Pendulum at t = 0 


Proposition 3.3 Reflections £ G permute the domains Mi and M 2 according to Tabled 


Proof Follows immediately from the definitions of the actions of reflections (2.201. 


□ 


Proposition 3.4 The domains Mi, M 2 are open, connected and simply connected. 


Proof From the definition of the sets Mi, M 2 (3.6l-(3.7l it follows that the domains Mj are homeomorphic to 

mS 


and therefore they are open, connected and simply connected. 


□ 


3.2 Decomposition in the Preimage of the Exponential Mapping 

We now consider the following set N G N corresponding to all potentially optimal geodesics: 

N = {{\,t)£N I t<t(A)}. 

By existence of the optimal geodesics, Exp(A'^) D M. In order to separate the Maxwell points in the preimage of 
the exponential mapping, introduce further the sets: 


N = nun', 

N'= {{X,t) £uI=iN, I 
TVj = NnN, i = 1,...,4. 


N = 




C Uj=iA^j 


t = t(A) or sin = o| U Ni, 


t<t{X), sin^ T^ojuiVs. 


3.2.1 Decomposition in N 


We now introduce the connected components Dj of the set N: 

N = DiU D 2 , 

Di = {(A,t)GuLiiV, I t<t(A), sin(^)>0}, 

D 2 = {(A,t)GuLiiV, I t<t(A), sin(^)<0}, 

where Di are defined explicitly in coordinates in Table(in the sets Ni, N 2 , N 3 ). Projections of the sets Di to 
the initial pha se cyl in der a re shown in Figure^ We note that for t < t (A ) = t ^^^(A) the values of p are given 
from formulas ( |2.29 (-(2.301, and the values of ti ^’^(A) are given in (2.221-( 2.241. The values of r in Table[^were 
calculated by using the definition of elliptic coordinates [T], formulas for Jacobi elliptic functions [?] and values 
of 7 and c from Figure[^ Note that enumeration of the sets Di is chosen to correspond to the sets Mj for further 
analysis. 

We now establish an important fact about the domains Dj that is vital in proving that the exponential 
mapping transforms Di diffeomorphically. 


Proposition 3.5 Reflections £ G permute the domains Di and D 2 as shown in Tabled 
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Di 

Di 

T >2 

A 

c? 




P 

(0, 2K) 

(0,2K) 

(0,2X) 

{0,2K) 

T 

{0,2K) 

{2K, AK) 

(2K, 4K) 

(0,2K) 

A 


C2 


C2 

P 

(0, 2K) 

iO,2K) 

(0,2K) 

{0,2K) 

T 

(0, 2K) 

i-2K, 0) 

(2K, 4K) 

{0,2K) 

A 

C^+ U 

ct u C 3 ^+ 

c^+uc^- 

0 

1*5 C 

1 

c 

0 

1- 

+ 

P 

( 0 , +00) 

( 0 , +00) 

( 0 , +00) 

( 0 , +00) 

T 

( 0 , +00) 

(—00, 0 ) 

(—00, 0 ) 

( 0 , +00) 


Table 2 Decomposition N = 


ld,e\e'=,e' 


Di 

D2 

D2 

Di 


Table 3 Action of on Dj C N 


Proof In paper we defined the action of reflections : N ^ N so that it satisfies the following properties: 


(A, t) = I o ^ 


ijf — jf 

ly — — ^ 




where ei i® pushforward of under the reflection eK Recall that eilly = —liy, for j = 1,2, 5, 6 

because these symmetries reverse the direction of time and e ity = ty, for j = 3,4, 7 because these symmetries 
preserve the direction of time [T]. Hence, it is sufficient to prove the case = D 2 as proof of all other cases 

{Di) = Dfc is similar. In order to prove the inclusion (Di) C D 2 we take any (A, t) = ( 7 , c, t) € Di and prove 
that 

: (A,t) H> (A^,t) = (7^,c^,t) G 02- 

By Proposition |3.1 1 

t(A^) = t o o (A) = t(A). 

Thus t < t(A). Moreover, at instant t/2 the trajectories of the vertical subsystem are given as: 

^t/2 = (7t/2!Ct/2) =e^‘'*^^(A), 


A ?/2 = 




_ J^yt /2 




Since A^ = o 


(A), we have 

,2 liyt/2 2 2 -liyt/2 2 ^ t / 2 / , \ 2/, N 

Aj /2 = e ' o e o e (Aj = e o e ' oe (A) = e o e ' (A) = e yXt/ 2 )- 


(3.8) 


In proof of (3.8 1 we used the fact that for any diffeomorphism F ■. M ^ M and a vector field on a manifold 
M, Fif^ = —F o o F. Clearly, £^{A(/ 2 ) = (7?/2>‘^t/2) from (6.3) [l] we have: 


(p't/2’^t/2'j — ( 7t/2Tt/2) • 


J, 

Thus sin= sin —< 0. We proved that (A^,t) C -D 2 , thus e^{D\) C D 2 . Similarly it follows that 


2 

£^( 7 ^ 2 ) C D\. Since o = Id, then o £^{Di) = Di => e^{D\) = 02 - 
Proposition 3.6 The domains Di,D 2 C N are open and conneeted. 


□ 


Proof Since —>■ is a diffeomorphism and {Di) = D 2 it suffices to prove that D\ is open and connected. 


Consider a vector field 


d 


d 


P = - ( c-^ — siny— ) € Vec(Af). 


97 


The flow of this vector field e^ is given as: 


= e^{\,t) = ( e 




= {'rt/2,C-t/2,t) ■ 
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Thus e^(-Di) = Di where 

5i = |(A,t)GiV I sin^>0, t<t(A)|. 

The set Di is a subgraph of a continuous function A i—>■ t(A) on an open connected 2-dimensional domain 
{( 7 , c) G C I 7 G (0, 27r), c G R}, thus Di is open and connected. Since Di = therefore Di is also 

open and connected. □ 


Proposition 3.7 There hold the inclusions: 

(1) Exp(A) C Mj, i=l,2, 

(2) Exp{N) C M, 

(3) Exp{N') C M'. 


Proof 


(1) It suffices to prove only that Exp(Di) C Mi, in view of the reflections eP Notice the decomposition: 

Di = (Hi n A^i) u (Hi n iV 2 ) u (Hi n N 3 ) u (Hi n N 5 ). ( 3 . 9 ) 

Let (A,t) G Hi n A^i = {(A,t) € Ni \ t < t(A), sin > 0} , thus p = I ^ {0,‘2K{k)). Further, from 
formula (5.3) [I] we have sisnr > 0. Now recall formula (3.2) [S]: 


sinh2t = Si 


2 k snp snr 

Z ’ 


A _ -1 7 2 2 2 

A = 1 — k sn psii r. 


(3.10) 


Then we get sinhai > 0, thus zt > 0, i.e., Exp(A, t) G M\. We proved that Exp(Hi n A^i) C Mi. All other 
required inclusions Exp(Hi ENj) C Mi, j = 2, 3, 5, are proved similarly, and the inclusion Exp(Hi) C Mi 
follows. 

(2) Since N = Hi U H 2 and M = Mi U M 2 , the inclusion Exp(N) C M follows from item (1). 

(3) We have N' = [N' n Ni) U (A^' n N 2 ) U {N' n N 3 ) U Na- 

Let (A, t) G N' D Ni = | (A, t) G Ni \ t = t(A) or sin = o| , then similarly to the proof of item (1) 

we get p = 2K{k) or snr = 0, thus 2 * = 0 by (3.101. From (3.6) [8] we get R 2 {qt) = dnr / 2 (p) ^ 0, 

and therefore ^ 0. We proved that Exp(A' n A^i) C M^ It follows similarly that Exp(A'' E Nj) C 

M', j = 2,3. Finally, if (A,t) G Na, then 


qt = {xt,yt,zt) = (t, 0,0) G M'. 


Consequently, Exp(A'') C M'. 


Theorem 3.1 For A G we have t“"'^(A) > t^^’^(A). 

Proof Apply equations ( |2.22[ )-( f2.24[ ) and ( |2.25[ )-( f2.28[ ). 
Proposition 3.8 The restriction Exp : N ^ M is non-degenerate. 


□ 


□ 


Proof From Theorem 3.1 tj°"^(A) > ti ‘^’'(A). Since for any v = [X,t) € N we have t < t(A) and therefore 


exponential mapping is non-degenerate 'du = {\,t) G N . 


□ 


Hence we proved properties PI, P2 and P3 for the exponential mapping Exp : H^ —> M^. It only remains to 
prove condition P4 now to establish that the exponential mapping Exp : Hj —>■ M^ is indeed a diffeomorphism. 


3.3 Diffeomorphic Properties of the Exponential Mapping 

In this subsection we prove that the exponential mapping Exp : H^ —>■ M^, i = 1,2, is proper. First we recall 
an equivalent formulation of the properness property. 

Definition 1 Let A be a topological space and {xn} C A a sequence. We write Xn —>■ dX if there is no compact 
K C X such that Xn & K for any n G N. 

Remark 1 Let A, Y be topological spaces and H : A —> T a continuous mapping. The mapping F is proper iff 
for any sequence {xn} C A there holds the implication: 

Xn^dX F{xn) dY. 

Below we apply this properness test to the mapping Exp : Hi —>■ Mi. 
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Lemma 1 Let {qn} C Mi. We have q-n —t dMi iff there is a subsequence {nfc} on which one of the conditions 
holds: 

(1) ^^0, 

(2) z —>■ +00, 

(3) a: —>■ 00 , 

(4) 2 / -t oo. 

Proof Any compact set in M\ is contained in a compact set € Mi | e < 2 < Ij/I < for 

some £ G (0,1). □ 

Lemma 2 Let {ir-n} C Di, then Vn — t dDi iff there is a subsequence {nff\ on which one of the following 
conditions hold: 


(1) 7t/2 

(2) 7t/2 27r, 

(3) Ct/2 oo, 

(4) t ^ 0, 

(5) t ->■ +00, 

(6) t(A)-t-^0 


Proof Any compact set in Di is contained in a compact set 

G iV I 7t/2 G [£,27r -e] , [cj/al < C [e, i], t(A) - t > e| , 

for some e G (0,1). □ 

Proposition 3.9 The mapping Exp : Di — t Mi, i = 1,2, is proper. 


Proof In view of the reflections , it suffices to consider the case Exp : Di —t Mi. Let {vn} C Di, Vn —t dD\, 
we have to show that qn = Exp(t'n) —t dM\. Taking into account decomposition (3.9l, we can consider the cases 
{vn} T D]_ C] Nj, j = l,2,3,5. 


Let {nn} C DinA^i, iy„ 
obtained in [T] and [5]: 


dDi. We will need the following formulas for the extremals At = e*^{A), A G Ci, 


'yf 

sin M = sifcsn((pf), 

^ = A:cn((pt), 

. fcsnpsnr . . , 2 2 2 

Sinn^:^ = si - — -, A = 1 — k sn psn r, 

R 2 {qt) = / 2 (p)Y^^dnr, / 2 (p) = dnpE(p) - k^snpcnp. 

Notice that p = |, 'r = V? + 5 , and consider all the cases (l)-( 6 ) of Lemmaj^ 

(1) If 7 t /2 Oj then sin 2^ = sifcsnr —>■ 0, thus sinhzf —^ 0, so zt —>■ 0, hence qn —>■ 9Mi(Lemmaj^ (I))- 

(2) If 7(^2 271, then sin 2^ = siksnr —>■ 0, thus sinh^t —>■ 0, so zt —>■ 0, hence qn —>■ dM\. 

(3) The case Ct /2 -^00 is impossible. 

(4) If t —>■ 0, then p —>■ 0, thus zt —>■ 0. 

(5) Let t —>■ + 00 , then p —> + 00 . Since p G (0, 2K(k)) then fc —>■ 1. Denote u = am(p) G (0, tt). On a subsequence 
we have u ^ u £ [ 0 i '^] and we will suppose so in the sequel. 

(a) If u G [0, tt), then p = F{u,k) ^ F{u, 1) = ; < + 00 , a contradiction. 

(b) Let u = ^, thus snp = sinw —>■ 1, cnp = cos(u) —>■ c. 

i. If snr —>■ 1, then Zi —>■ 0, thus zt —^ 00 . 

ii. Let snr —>■ s 1, then dnr —>■ \/l — yl: 0. Denote 


<? 2 (w) = f 2 {F{u, k)) = \/1 — k'^ sin^ uE(m, k) — k^ sin(u) cos(u). 

We prove now that —>■ + 00 , then —>■ + 00 , thus R 2 {qt) °o, so +yi + zf —>- 00 , whence 

qt —>■ dMi. Denote k' = %/! — —>■ 0. We can suppose that on a subsequence —>■ a G [0, + 00 ]. 

We have 


k'^ sin(u) cos(u) = sin(u) cos(m) + o{k'^), 

Vl - sin^ u = \J COS^ U + — A;'2 


COS^ U. 














12 


Yasir Awais Butt, Yuri L. Sachkov, Aamer Iqbal Bhatti 


Now we estimate E{u, k) from below: 

u u u 

E(u,k) — s\n{u) = J \/1 — k'^ sin^ tdt — J cos{t)dt = j 


1 — sin^ t — cos^ t 
\/\ — k'^ sin^ t + cos t 


dtt 


1 - k^ 


1 - k^ 


sin^ t dt 


sin(2M) 


= ^fc'2(l + o(l)). 


Thus, 


E{u, k) > sin(t6) + + o(l)). 


A. Let a € [0, +oo). Then cos(m) = ak' + o(k'), sin(u) = 1 + o(l), thus 

k^ sin(u) cos(«) = ak' + o{k'), 

^ 1 — k'^ sin^(M) = \/l + a'^k' + o{k'), 

E{u,k) = l + o(l), 

y/1 — k^ sin^ uE(u, k) = \/1 + a^k' + o{k'), 

92 {u) = (^\/l + k' + o{k'), 


92 


(u) _ (Vl + a^ - a) 


fc '2 

and the claim follows. 

B. Let a = +oo, thus k' = o(cos(m)). Then 


k' 


(1 + o(l)) “t oo, 


k^ sin(M) cos(u) = sin(u) cos(u) — k'^ cos(u) + o cos(u)^ 

\/l — u = cos(u) Win- ^ -h o f ^ 

y cos2 u y cos21 


1 

= cos(m) + --^ + o 


2 cos(u) ycos(t6) 

/-^- 1 f 

V 1 — sin^ uE(u, k) > cos(m) sin(u) + --+ o ( 

2 cos(m) \ 


cos(m) 


92{u) > ^^(l + o(l)), 

+ °(1)) +°0’ 

and the claim follows. 

hi. Let u G (0,7r), then / 2 (p) = 92 {u) —t Icosttl {E{u, 1) + sinu) > 0, thus 


f2{p) 


+ 00 . 


Since > 1, then i? 2 (?t) —t oo, so a;f + t/f + 2 ^ —t oo, whence qt —t dMi. 

iv. If u = TT, then snp = sin(ti) —>■ 0, thus zt —t 0. 

(6) Let t(A) — t —>■ 0. Recall that t(A) = AK(k) for A G Ci, thus AK(k) — t —>■ 0. Since k G (0, 1), then there is 
a subsequence {rim} on which k ^ k £ [0,1]. If fc G [0,1), then K{k) —^ K{k) < +oo, thus t —>■ AK{k), so 
p = 2K{k). Consequently, sinh 2 t —>■ 0, whence qn —t dMi (Lemma[^ (I))- If fc = 1, then K{k) —^ +oo, thus 
t —> + 00 , qn —t dMi by Item (5). 

Consequently, in each of the cases (l)-(6) of Lemma we get —>■ dMi for a sequence {vn} C Di n 

A^i, Vn —r dDi. All the rest cases {vn} C Di n Nj, j = 2, 3, 5, are considered similarly. 

Summing up, for any sequence {vn} C D\ with Vn —>■ dD\ we have Exp(i/n) —> dMi. Thus the mapping 
Exp : Di —>■ Ml is proper. □ 

Now we get the main result of this section. 
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Theorem 3.2 The mapping Exp : Di —> Mi, i = 1,2, is a diffeomorphism. 
Proof All of the conditions PI—P4 are satisfied for the mapping Exp : Di —> Mi: 


— Di C N and Mi C M are open subsets thus 3-dimensional manifolds (Proposition 3.6 Proposition 3.41, 

— PI - Di is connected (Proposition |3.6| , 

— P2 - Ml is connected and simply connected (Proposition |3.4| |, 

— P3 - Expln^ is non-degenerate (Proposition |3.8 [) , 

— P4 - Exp : Di —> Ml is proper (Proposition |3.9p . 

Thus Exp : Di —> Mi is a diffeomorphism. By virtue of the reflections, Exp : D 2 —t M 2 is a diffeomorphism as 
well. 

□ 

Corollary 1 The exponential mapping Exp : N ^ M is a diffeomorphism. 

Proof Follows from Theorem |3.2[ □ 


3.4 Cut Time 


Now we can prove that inequality ( |3.5[ ) is in fact an equality for A G C\C' 4 . 
Theorem 3.3 //A € C\C' 4 , then tcut(A) = t(A). 


Proof Let A € C'\C '4 = uf_iCi U C 5 . In view of inequality ( 3.5 I, it remains to prove that tcut(A) > t(A). Take 
any ti € (0,t(A)).We need to prove that the geodesic Exp(A,t) is optimal on the segment t G [0,ti]. 

Consider first the case A G uf—iCi. If sin '''‘2^^ ^ 0, then (A,ti) G N, and gi = Exp(A,ti) G M. By virtue of 
Proposition 3.7 and Theorem 3.2 the point gi has a unique preimage under the mapping Exp : N M. Thus 


the geodesic Exp(A,t) is optimal on the segment t G [0,ti]. 

If A G uf—iCi and sin '^* 2 ''^ ~ then we can choose t 2 G (ti,t(A)) such that sin 7 ^= 0. By the 

argument of the preceding paragraph, the geodesic Exp(A, t) is optimal at the segment [0, 12 ], thus at the segment 
[0, ti] C [0, 12 ] as well. 

Finally, if A G C 5 , then (A,ti) G N, and the geodesic Exp(A,t), t G [0,ti], is optimal as above. 

We proved that tcut(A) > t(A), thus tcut(A) = t(A) for any A G C'\C 4 . □ 

We will be able to prove the equality tcut(A) = t(A) for X G C 4 below after the description of the structure 
of the exponential mapping Exp : N' —>■ M^ The geodesic Exp(A,t), A G ( 74 , requires a separate study since it 
belongs to the set M' for all t > 0. 

establishes the fact that since Exp : A( —>■ M is a diffeomorphism, hence upto time 


Intuitively, Theorem 


3.3 


t < t(A) there is a unique point v = {X, s) £ N that is mapped to a unique extremal trajectory qs = Exp(A, s) G 
M that joins qo £ M to qi £ M d M. Hence, the trajectory qs = Exp(A, s) £ M is optimal and therefore 
tcut(A) = t(A). It therefore follows that optimal synthesis in the domain M is given by: 

Uiiq) = hi{X), i = l,2, (A, t) = Exp“\g) G fV, g G M, 

where Ui are the control variables (i.e., translational and rotational velocities) and hi are the optimal controls 
defined in (4.8) [I]. 


4 Exponential Mapping on the Boundary of Diffeomorphic Domains 

Until now we have studied the mapping Exp : N ^ M and proved that it is a diffeomorphism. This allowed us 
to prove that the cut time tcut(A) = t^'^’'(A), A G C\Ca. In this section we obtain the global structure of the 
exponential mapping in order to characterize the cut locus and the Maxwell strata and to construct the optimal 
synthesis. Specifically we study the mapping Exp : N' —>■ M' where: 

A(' = {(A,t) G I t = tf^^{X) or sin(^) =0}u{(A,t) G Ar 4 | t < 27r = t“"j(A)} , 

M' = |gGM I 2 :^-by Vo, « = o| . 
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3 

A 

p 

T 

k 

1 

Cl* 

2K 

(0,Y) 

(0,1) 

9 

ct 

2K 

(0,Y) 

(0,1) 

17 

cv 

2K 

K 

(0,1) 

21 

cy 

2K 

0 

(0,1) 

25 

Ct 

2K 

0 

(0,1) 

29 

Ct 

2K 

K 

(0,1) 


Table 4 Decomposition Yj, j S {1, 9,17, 21, 25, 29} 


A 

P 

T 

k 

cy 

(0,2Y) 

0 

(0,1) 


(0,2Y) 

0 

(0,1) 

_ 

(0, + 00 ) 

0 

1 


Table 5 Decomposition Yj, j = 35 


3 

A 

t 

33 

cy 

27 r 

39 


(0, 27 r) 


Table 6 Decomposition Yj, } £ {33,39} 


C 



Fig. 4 The sets Yj with t = or sin = 0 


4.1 Stratification of N' 


We define subsets NjCN', ji = 1,..., 40, as follows: 

— for j £ {1, 9,17, 21, 25, 29} the sets are given by Tabled for j = 35 by Tablelsland for j £ {33, 39} by 
Table § 

— for all the rest j the set Y, are defined by the action of reflections e* as in (4.1 l“(4.4l: 



II 

+ 

i = 1,... ,7, j = 1,9, 

(4.1) 

(Mr) 

= 

i = 1,2,3, 

(4.2) 

(Yj) 

II 

+ 

i = l,2,3, j = 21,25,29,35, 

(4.3) 


= N'+i, 

j = 33,39. 

(4.4) 


The following stratification of the set N' follows from the definition of the sets Yj. 

Lemma 3 The stratification of N' shown in Figures I4I5I is given as: 

Y' = uj°iY'. (4.5) 


From Figures 4 5 we see the sets Y) given in Tables Ul 5l 6 pertain to the quadrant of the phase portrait of 


vertical subsystem for which A = ( 7 , c) £ C such that 7 ^llTw] and c £ [0 , 00 ). For A = ( 7 , c) in other parts of 
phase portrait, the sets Y) are obtained by the reflection symmetries (4.1l-(4.4l of the vertical subsystem. 










































Cut Locus and Optimal Synthesis in Sub-Riemannian Problem on the Lie Group SH(2) 


15 



Fig. 5 The sets N’- with t < if sin ^ = 0 



Fig. 6 Stratification of the quadrant Q 


4.2 Stratification of a Quadrant of the Plane z = Q 

Define the following curves and points in the quadrant Q = {(x, y) G | x > 0, t/ < O} (see Figure j^: 


71 : 

72 : 

73 : 

74 : 

75 : 

P : 

O : 


X = 0, y = yi{k) = - 


4a(fc) 
Vl - 


, fee (0,1), 


,,. Aka{k) , 4a(fc) , m 

X = X2(k) = ^_\.2 ’ y = = ~ i _fc2 ’ fee (0,1), 

4 Ak 

X = X3(fc) = j-^P(fc), y = y3{k) = -^-^E{k), fee (0,1), 

X = X 4 (f) = t, y = 0, t G (0, 27r), 

4 


X = 27r, y = 0, 
X = 0 , y = 0, 


E{k), y = 0, fcG(0,l), 


where a(k) = E(k) — (1 — k^)K{k), k G (0,1). The curves 71 ,... ,75 result from substitution of t = f^“(A), 
and ip = T — p from Table in the equations of extremal trajectories for A G The curves 71 ,..., 75 and 

the point P are the images of certain sets Exp {Nj'j under the projection 


p:{q£ M \ 


2 = 0 }-s-(x,i/, 0 ) i-i-(x, 2 /). 


(4.6) 
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71 =po Exp (A^ag) . 

72 =p o Exp {N 25 ) , 

73 =po Exp {N 21 ) , 

74 =P ° Exp (A^ag) , 

75 =p o Exp (A^^y) , 

P =po Exp (A^aa) . 

These equalities can be verified easily. From |8j we know that the first Maxwell points with t = and 

conjugate points with t = (A) and snr cnr = 0 lie in the plane 2 = 0. Hence, the curves 71,..., 75 decompose 

the fourth quadrant of the plane 2 = 0 into various regions (see Figure]^. The regularity and mutual disposition 
of the curves 71,..., 75 are described in the following lemmas. 

Lemma 4 The function a{k) satisfies the following properties: 

a : (0,1) —>■ (0,1) is a diffeomorphism, (4.7) 

A: —>■ 0 =7 a(fe) = + o{k^), (4.8) 

k^l -0 =7 a{k) = 1 - In + 0(fc'2) (4.9) 

where k' = %/! — Moreover, the function a(k) is convex. 


Proof If fe —>■ 0, then 


E{k) = lU-’^]+o{k^), 


which gives asymptotics (|4.8[|. If fc —>■ 1 — 0, then 


K{k) = In +o{k'), 

E{k) = l + +0(fc'"), 


which gives asymptotics (4.9l. Finally, property (4.71 follows since 


| = /,RW>0, 

lim a(fe) = 0, 
fc—>-0 

lim a(k) = 1. 
fc^i-o 

The function a(fc) is convex since ^ = kK{k) increases Vfc £ (0,1). 

Lemma 5 The function y = yi{k) defines a diffeomorphism yi : (0,1) —^ (—oo,0). Moreover, 

lim yi{k) = 0, 
fc—>-0+ 

lim yi(k) = — 00 . 
fc-n- 

Proof The function y = yi{k) is a strictly decreasing function with: 

dyi _ -AkE{k) 


□ 

(4.10) 

(4.11) 


dk (I_jfc2)5 

Further, Lemmayields the asymptotics: 


<0, fee (0,1). 




1-0 

and the statement of this lemma follows. 


Vi(.k) -pr -00, 


□ 
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Lemma 6 The function x = X 4 {t) defines a dijfeomorphism X 4 : (0, 27r) —>■ (0,27r). Moreover, 

lim xMt) = 0, 

t-i.0+ 

lim X 4 ,(t) = 2 'k. 
k—^2'7z~ 

Proof Clearly X 4 {t) is a smooth bijection with a smooth inverse. Hence it is a diffeomorphsim. The limits can be 
calculated by direct substitution in X 4 {t). □ 

Lemma 7 The function x = x^{k) defines a diffeomorphism X 5 : (0,1) —t (27r, + 00 ). Moreover, 

lim xsik) = 27r, 
fc—> 0 + 

lim X 5 {k) = + 00 . 
fc-s-i- 

Proof The function x = x^{k) is a strictly decreasing function with: 

dx 5 ^ 4a(fc) ^ p 
dk k{l-k‘^)l 


k ^0 E{k) 


xz{k) —>• 2-k, 


fc —t 1 — 0 E{k) — 1 - 1 xz{k) —>• + 00 , 


and the statement of the lemma follows. □ 

Lemma 8 The functions x = X 2 {k), y = y 2 {k) k € (0,1), define parametrically a function x = X 2 {y) which 
is a diffeomorphism X 2 ■ (— cx3,0) —t (0, + 00 ) with limy_>._oo * 2 ( 1 /) = + 00 , limj,_j.g- X 2 {y) = 0. Moreover, 

-y-2 <x 2 {y) < -y, y € {-oo, 0 ). (4.12) 

The curve 72 is convex, has near the origin the asymptotics 

y =—Tv^x^ + o , X —i 0, (4.13) 

and has an asymptote y -\- x -\- 2 = Q as x ^ 00 . 


Proof Notice that 


fc —t 0 => X2{k) —t 0, y2{k) —0, 

fc — t 1 => X 2 {k) — t + 00 , J/2(fc) —t —CX3. 


dx 2 _ 4 ((1 + fc") E{k) - (1 - e)K{k)) _ 4 (a(fc) + k^E{k)) ^ 

dk (l-fc2)2 k{l-k'^f 

dy^ _ Ak{2E(k)-{l-k^)Kik)) _ 4fc(a(fc) + L(fc)) 

dk (l-fc2)2 { 1 -k'^f 

thus the functions X 2 {k) and y 2 {k) define diffeomorphisms X 2 '■ (0,1) —>■ (0, + 00 ) and j /2 : (0,1) —>■ {— 00 , 0). So 
these functions define parametrically the diffeomorphism 

x = X 2 {y), j/G(-oo,0 ), a;€(0,+ 00 ), 

y = y 2 (,x), ®g(0 ,+ 00 ), ye (- 00 , 0 ). 


lim X 2 {y) = lim X 2 {k) = + 00 , 

y—>■—00 fc —>1 

lim X 2 {y) = lim X 2 (k) =0. 
y-> 0 - fc-> 0 - 


Notice that 
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Fig. 7 The curve 72 and its bounds y -\- x = —2, y + x = 0. 


Now we show that the curve 72 is convex. We have 
dj/2 _ dy2/dk 


dx dx 2 /dk 


= a{k), 


a{k) = —k 


da 

dk 


2E{k) - (1 - k^)K{k) 

(1 + k‘^)E{k)- (1 - k^)K(k)’ 

2 \ 3E'^{k) - (5 - k^)E{k) K{k) + 2(1 - k^)K‘^{k) 




((1 + fe 2 ) E{k) - (1 - A;2) K{k)Y 


(4.14) 

(4.15) 


Since a(k) = E{k) — (l — K(k) G (0,1), then ^ ((^ “ numerator of the function 

t = I—>■ — (5 — t + 2 (l — is negative for t £ ((l — , l) thus the numerator of fraction (4.15 I 

is positive. Therefore, ^ > 0, i.e., is increasing for k £ (0,1) and also increasing for x £ (0,+oo). Thus the 
function y 2 {x) and its graph, i.e., the curve 72 , are convex. The second inequality in (4.121 follows since 

X2(k) 


1/2 (fe) 


= -fe>-l, fcG( 0 ,l). 


The first inequality in (4.12 1 and existence of the asymptote y + x + 2 = 0 follows from equalities: 

r vYk) 1 
,lim —jTX = - 1 , 
fe->l- X 2 {k) 

lim {y 2 {x) + X 2 {y)) =- 2 , 

fe— ^1 — 

( 2 / 2 ( 2 :) + 2 : 2 ( 2 /)) + 2 = (1 + fe - 2 a(fc)) > 0 , 


since a{k) < k < for k £ (0,1). Finally asymptotics (4.131 follows since 


X 2 {k) = Tvk^ + o{k^), y 2 (k) =—nk'^ + o{k‘^), fc —>• 0 . 


□ 

A plot of the curve 72 with its bounds given by ( |4.12[ ) is shown in Figure]^ 

Lemma 9 The functions x = ® 3 (fe), y = y^{k), define parametrieally a function x = 2 : 3 ( 2 /) which is a diffeo- 
morphism 2:3 : (—oo,0) —/ (27r, + 00 ) with limy_)._oo 2 : 3 ( 2 /) = + 00 , limj^_^Q+ 2 : 3 ( 2 /) = 27r. Moreover, 

2 : 3 ( 2 /) >27r, 2 : 3 ( 2 /) > 2 - 2 /, 2 /e (- 00 , 0 ). (4.16) 


The curve 73 is convex and has an asymptote y + x = 2 as x ^ 00 . 
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Fig. 8 The curve 73 and its bounds y + x = 2 , x = 27 r. 


Proof Notice that 

k —>■ 0 X 3 (k) —t 27r, y 3 (k) —t 0 , 

k —>■ 1 a: 3 (fc) + 00 , ysik) —t — 00 . 

Furthermore, 

dx3 _ 4 ((1 + fc") E{k) - (1 - k^)K{k)) _ 4 (fl(fc) + k^Ejk)) ^ 
dk k{l-k'^f k{l-k'^f ^ ’ 

dm ^ i{2E{k)-{l-k^)K{k)) ^ 4(a{k)+E{k)) ^ ^ 

dk k{l-k^f k(l-k‘^f 

thus the functions ® 3 {fe) and y 3 {k) define diffeomorphisms *3 : (0,1) —t {27r, + 00 ) and 1/3 : (0,1) —^ (— 00 ,0). So 
these functions define parametrically a diffeomorphism 

X = X 3 {y), t/G(-oo,0), a: G (27r,+ 00 ). 


Notice that 


lim X 3 {y) = lim X 3 {k) = + 00 , 

y^—00 fe— >1 

lim X 3 (v) = lim xsik) = 2tt. 
v^0+ ' ' k^0+ 


Since > Q, therefore X 3 {k) > 27r for k G (0,1), which gives the first inequality in (4.16 1 . The second inequality 
in (4.161 and existence of the asymptote y + x = 2 follow from the equalities: 

lim ^ = -1 

X3(k) 

lim (1/3(2;) + 2:3(1/)) = 2 , 

fc— >-1 

(1/3(2:)+ 2:3(1/))-2 = 


1 + fc 


E{k)- 


1 + k 


> 0 . 


Finally, convexity of the curve 73 follows since 

dm ^ dy3/dk 
dx dx 3 /dk 


= aik), 


where a{k) is given by (4.141, which is increasing by the proof of Lemmal^ 
A plot of the curve 73 with its bounds given by (|4.16[) is shown in Fig^ 


□ 


Lemma 10 For any y G (—oo, 0 ), we have 2:2(1/) < 2:3(1/). 
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j 

y 

X 

z 

1 

(—00, 0) 

(x 3 (y),+oc) 

0 

9 

(—00, 0) 

(0,x2(y)) 

0 

17 

0 

( 277 , +00) 

0 

21 

( —CXD, 0) 

X 3 {y) 

0 

25 

(—00, 0) 

X2(y) 

0 

29 

(—00, 0) 

0 

0 

33 

0 

27r 

0 

35 

(—00, 0) 

{x 2 iy),X 3 {y)) 

0 

39 

0 

(0,27r) 

0 


Table 7 Definition of C p ^(Q). 


Proof It follows from Lemmas]^ and l^that X 2 {y) < -y < 2 - y < X 3 {y), y G (-oo, 0). 
Lemmas [s ■ 10 allow us to define the following domains in the plane Q C . 


^l,y- 


mi = j 


1 y <0, 

0 <x < X2{y)'^ , 

m2 = \ 

^^{x,y) G 

1 y <0, 

X2{y) <x < a; 3 (t/)|, 

m3 = \ 

^^{x,y) G 

1 y <0, 

X 3 {y) < , 


□ 


see Figure 

Lemma 11 The domains mi,m 2 ,rn ,3 C ^x,y 0,'’'^ open, conneeted and simply connected, with the following 
boundaries: 


dmi = 7 i U 72 U {O}, 
dm 2 = 72 U 73 U 74 U {O, P}, 
dm 3 = 73 U 75 U {P}. 

Moreover, the quadrant Q has the following decomposition into disjoint subsets: 

Q = U (uf=i 7 ,) U {0,P}. 

Proof Follows from the definition of the domains m^ and from Lemmas [sHlOl CH 


Define the inverse images of the sets mi,yi, and P via the projection p (|4.6[): 


Mg = p“^(mi). 

M33 =p ^(m2). 

M'i =p ^(m3). 

M29 =p~^('yi), 

M25 =P“^(72), 

M21 =P~^{T 3 ), 

M39 =p"\74), 

Mir=p-\-/ 5 ), 

M'33=p-\P). 

Explicitly, these sets are defined in Table 




Now we aim to prove that all the mappings Exp : Nj —>■ Mj are diffeomorphisms for the sets Nj and 
dehned by Tables 0ii[Z] 

Lemma 12 For any j G {17, 21, 25, 29, 33, 39} the mapping Exp : Nj —>■ is a diffeomorphism. 

Proof Follows immediately from above lemmas: 

— Lemma 

— Lemma 

— Lemma 

— Lemma 

— Lemma 

— and it is obvious for j = 33 . □ 

Now we consider the mappings of 2-dimensional domains. 

Lemma 13 The mapping Exp : Ng —>■ Mg is a diffeomorphism. 


for j = 17, 
for j = 21, 
for j = 25, 
for j = 29, 
for j = 39, 
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Proof In the coordinates p = ^ and r = (ip + §) /^, the domain Ng is given as follows: 

iVgiAGC^, S2=0, p = 2K{k), TG{0,K{k)), fee (0,1). 

Introduce further the coordinate u = am(T), then, 

Ar':s 2 = 0 , p = 2K{k), u G (o, , fcG(0,l). 

In these coordinates the exponential mapping Exp(A, t) = (x, y, z) is given as follows: 

, 4fca(A:) cos(u) 

X = xg{u,k) = - i_f.2 -’ 

, 4a(fc) \/l — sin^ (u) 

y = yg{u,k) = -p-, 

2 = 0 . 

Consider the mapping: 

/q ■ ^u,k ^ (^ 95 ^ 9 )? 

Du,k = (o,|)^ X (0,l)fe. 

We have to show that the mapping /g : D —t mi is a diffeomorphism. 

(1) First we show that fg{D) C mi. 

We fix any k G (0, 1) and show that the curve F : u —>■ (xg,yg), u G (O, ^) , is contained in mi. Compute 
hrst the boundary points of P: 


Further, since 


u ->• 0 r{u) ->• {X2{k),y2{k)) G 72 , 
u -i- ^ =7 r{u) -)■ ( 0 ,i/ 2 {fe)) G 71 - 


dxg 

du 

dyg 

du 


Aka(k) 

1 - fe2 


sin(u) < 0, 


4:k'^a(k) sin(M) cos(m) 

^ sin^(u) 


then the curve F is a graph of the smooth function x 1 —>■ yg{x). Since 

dygfdu _ kcos{u) 


dyg ^ _ 

dx dxg/du 

then the curve F is concave. Moreover, 


•\/l — sin^(M) 


for u G 



dyg 

dx 


u=0 


—k > a{k) 


dm 

dx ’ 


where a{k) is given by (4.141. Since the curve 72 is convex, it follows that the curve F lies below the curve 
72 . Thus F C mi. Consequently, /g(F) C mi. 

(2) Since 

d{xg,yg) _ l 6 k^E{k)a{k)sin{u) 


d{u,k) (^i _ 

sin^('u) 


>0, 


(4.17) 


then the mapping /g : F —>■ mi is non-degenerate. 

(3) Finally we show that the mapping /g : F —>■ mi is proper. 

It is obvious that a sequence {un, kn) ^ dD iff it has a subsequence on which at least one of the conditions 
hold: 

M —>■ 0, u—A: —> 0, A: —>• 1. (4.18) 


On the other hand, a sequence {xn,yn) —t dm\ iff it has a subsequence on which at least one of the conditions 
hold: 


*—>■0, a; —>• -boo, y—>■ 0, j/— >■ — 00 , X 2 (y) — x ^ 0. (4.19) 

We show that in each of the cases ( |4.18| l we have one of the cases ( |4.19| |. If A: —>■ 0, then xg —>■ 0 and t/g —>■ 0. 
We can assume below that A: —>■ A G (0,1]. 

Let fcG(0,l).IfM—>0, then (xg,yg) —>■ (x 2 (k), y 2 (k)) G 72 thus X 2 (y) — x —>■ 0. If u —>■ then a;g —>■ 0. Let 
fc = 1. If u —>■ 0, then xg — >■ 00 . If u —>■ then yg —>■ 00 . 

We proved that the mapping /g : F —>■ mi is proper. 
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(4) The sets D, mi C are open, connected and simply connected. 

Thus fg:D^ mi is a diffeomorphism, as well as Exp : Ng —>■ Mg. 

Lemma 14 The mapping Exp : A'’( —>■ is a diffeomorphism. 

Proof In the coordinates p = § and r = p + |, the domain is given as follows: 

: A G C?, Si = 0, p = 2K{k), r G (0, K{k)), k G (0,1). 

Introduce further the coordinate u = am(T), then 

N[-.si=0, p = 2Kik), ug( 0,|], fcG(0,l). 

In these coordinates the exponential mapping Exp(A, t) = [x, y, z) is given as follows: 


X = xi(u, k) 

y = yi{u, k) 
2 = 0. 


AE{k) Y^l — sin^(u) 

1 - fc2 

AkE(k) cos(m) 


Consider the mapping: 


/l ■ ^u,k ^ (^7^) ^ 

Du,k = (o.0^ X (0,l)fe. 

We have to show that the mapping /i : D —t m 3 is a diffeomorphism. 

(1) First we show that fi{D) C m 3 . 

If {u,k) G D, then xi{u,k) > 0, yi{u,k) < 0, thus fi{D) C = {(a:,y) G | a: > 0, 

The boundary of the domain m 3 in R^_ is the curve 73 and along this curve we have = —k. 


73 = < {x, y) G Re¬ 


consider the function 


m3 = { {x,y) G R+_ 


X > 


1-4 

1-4 


(fii {u, k) = X — 


1 - 


x=xi{u.,k), y=yi {u,k) 


We have to show that ipi{u, k) > 0 for (u, k) G D. Since 

AE{k)^/l - sin^iu) AE{k) 


y>i(u,k) = 


1 - 

4\/l — sin^(ai) 

1 - A:2 


1 - 


k^ cos^ u 
1 — k^ sin^ u 


E(k) - E{k)ffl - fe2 sin^(M) 


where k = have to show that 

V 1 —sin^ u 


V>2(u,k) =E{k)-E{k) 



fc2 sin^{u) > 0, 


(u, k) G D. 


Since P 2 ( 0 , k) = 0 and 


d(p2 ^ 
du 

where <fi 3 {u,k) = (l — k^ sin^{u)'j E(k) — (l 
{u, k) G D. By Lemmawe have 


tan(u) 


zlfi3{u,k), 


\/l — A ;2 sin^(M) 

— k^\K{k), it is sufficient to show that 7 ) 3 ( 14 , fc) > 


y < 0}- 

Thus 


0 for all 


a{k) = E{k) - (l-k'^'^ K{k) > 0, fcG(0,1), 
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thus 

a(k) = E{k) - (l - K{k) 

(1 - sin2(u)) E{k) - (l - k^) K{k) 
1 — sin^(u) 


That is, ip 3 {u,k) > 0, 
(2) Since 


V(m, fe) G D. Thus it follows that fi{D) C m 3 , i.e., Exp(Af() C M(. 


d(xi,yi) 
d(u, k) 


\&E{k) a{k) sin(u) 

(1 


< 0 , 


then the mapping /i : D —t m 3 is non-degenerate. 

(3) Finally we show that the mapping /i : D —t m 3 is proper. 

In order to show that the mapping /i : D —^ m 3 is proper, we show that if a sequence (un, kn) £ D satisfies 
one of the conditions: 

TT 

u —t 0 , u —t , k —^ 0 , k —t 1 , 


then its image {xn,yn) = fi{un,kn) satisfies one of the conditions: 


2 ; —> 0 , ® ^-boo, y —t 0 , y —t 00 , X 3 {y) — X ^ 0. 


We can assume that fc —>■ fc G (0,1], m G u G [0, ^]. If fc = 0, then yi —>■ 0. 

Let k G (0,1). If It —>■ 0, then ( 2 : 1 , yi) —>• {x 3 {k),y 3 {k)) G 73 , thus 2 : 3 ( 1 /) — 2 ; —>• 0. If ft = 5 , then yi —>■ 0. Let 
fc = 1. If u G [0, 5), then a;i —>■ 00, yi —>■ 00. Let u = ^, then 


yi 


XI 


4cos(u) 

■ l-F ’ 


1 


1 - 




COs( 2 t) 
1 - 


We can assume that —>■ d G [0, -boo). If d G [0, + 00 ), then 2:1 —>■ -boo, and if d = -boo, then yi —>■ 00 . 

We proved that the mapping /i : D —>■ m 3 is proper. 

(4) The sets D, m 3 C are open, connected and simply connected. 

Thus f\ ■. D ^ m 3 is a diffeomorphism, as well as the mapping Exp : A'’( —>■ M[. □ 

Lemma 15 The mapping Exp : N 33 —> M 35 is a diffeomorphism. 

Proof It follows from Tables that 

N 35 = {i\t)(^N I 7|=0, C|>0, tG(0,t(A))}, 

M 3 ^ = {q£M I 2 = 0, y < 0, 2 : 2 (y) < 2 ; < 2 ; 3 (y)} . 

Further we have an obvious decomposition 

N35 = A^ss,! LI N 33^2 U N33^3, 

=^35niV,, i = 1,2,3. 

(1) We show hrst that Exp(A^ 35 ) C M 33 . 

Consider the set N 3 ^ 2 - In th® coordinates p = k ^nd r = (p + 7 ) /k, the domain A^ 35^2 is given as follows: 

Af 35.2 : A G C 2 +, S 2 = l, p={Q,2K{k)), r = 0, A: G (0,1). 

Introduce further the coordinate u = am(p), then 

A^ 35,2 : A G C 2 +, S 2 = l, u=(0,27r), r = 0, fc G (0,1). 

In these coordinates the exponential mapping Exp(A,t) = (x,y,z), (A, t) G 2 is given as follows: 


sin(u)y 1 — sin^( 2 t) — cos(u) a(u, k) , 

y = y35{u,k) = 

2 = 0 , 


\Jl — k"^ sin^ (u) ot{u, k) — k^ sin(M) cos(u) 


^ ^ 2k 

X = X35(u,k) = 
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where a{u,k) = E{u,k) — (l — F{u,k). Thus Exp{Af 35 2 ) C {q € M | a = 0}. Now we show that 

X 35 (m, k) > 0, y 35 {u, k) < 0 for (u, k) € (0, J) x (0,1). We have to prove the double inequality 


k) < a{u, k) < Q 2 (u, fc), (u, k) G ( 0 , x ( 0 , 1 ), 


ai{u,k) 

a2{u,k) 


k^ sin(M) cos(u) 

•\/l — sin^(M) 

sin{u) sin^(u) 

cos(u) 


This double inequality follows since 


Qi(0, k) = a(0, k) = 02 ( 0 , k) = 0, 
^ {a{u, k) — a\(u, k)) = — k^^ sin^(u) > 0 , 

^ (q2(w, k) — a{u, k)) = 1 — k^ > 0. 


Thus X 35 {u,k) > 0, y 35 {u,k) < 0 for (u,k) G (O, 5) x (0,1). If u G [5,7r), k G (0,1), then sin(u) > 

0, cos(ti) < 0, a{u,k) > 0, thus X 35 {u,k) > 0, t/35(w, fc) < 0. We proved that 

Exp(A ^35 2 ) C {q G M I 2 = 0, a; > 0, y < 0}. The sets N 35 i and 3 are considered similarly. 
Thus it follows that 


Exp(Ai"35) C _:= {g G M | 2 = 0 , ® > 0 , 1/ < 0 } . 


We now show that Exp(A^ 35 ) C Notice the decomposition 

M+_ = U Mg U M 21 U M 25 U M 35 . 


By contradiction, let Exp(Al 35 ) M 35 , then Exp(N 35 ) n (M{ U Mg U M 21 U M 25 ) 9 q. Let q G Exp(N 35 ) n 
M[ (the cases of intersection with Mg, M21, M25 are considered similarly). Then there exist (A35, tss) G N35, 
(Ai,ti) G N[ such that q = Exp(A35,t35) = Exp(Ai,ti). Notice that 


(A35,l35) e IV 35 
(Ai,ti) G Ni = 


^ l35 < Icut (A 35 ) , 
tl < tcut (Al) . 


(4.20) 

(4.21) 


If 135 < fi, then the trajectory Exp (Ai, t), t € [0, t\], is not optimal which contradicts to (4.211 . If 435 > ti, 
then the trajectory Exp(A 35 ,t), t G [ 0,135 + e] is not optimal for small e > 0 which contradicts to (4.201. 
Thus Exp(A^ 35 ) n M( = 0. Then it follows that Exp(N 35 ) C M 35 . 

(2) We now prove that Exp : N 35 —>■ M 35 is non-degenerate. 

Let V = (A,l) G N' 2 ,^ 2 - In fhs coordinates (jp,T,k) on N 35 2 , we have p G (0,2iL(fc)), r = 0, fc G (0,1). 
Since t < 4if(fc) = lcut(A) < 1^°’^^(A), therefore the Jacobian ^{v) 0. We have 


dq_ ^ d{x,y,z) 
du d{p, r, fc) 


Xp Xt X]^ 

yp yr yk 

Zp Zr 2fe 


Since Exp (lVi_2) C {g G M | 2 = 0}, then Zp{u) = Zj~{u) = 0, thus 


dq 

dv 




d{x,y) 
d{p, fc) 


{v) Zr{v) ^ 0 , 


so {v) ^ 0. Since v G A^35,2 is arbitrary, then Expiry/ is non-degenerate. Similarly it follows that 

Exp is non-degenerate at any point v G A ^35 1 U A ^35 3 . 

(3) The mapping Exp : AI 35 —^ M 35 is proper. This follows similarly to the proof of properness of Exp : D\ —> Mi. 

(4) It is obvious that M 35 is a connected, simply connected 2-dimensional manifold. In order to prove the same 
property for AI 35 , consider the vector field 


7 ^ = Cif— siny^ G Vec(Al). 
07 ac 


Since 

e*/ 2 ^(Ar^ 5 ) ={(A,l)GAr I 7 = 0 , c>0, 1 < t(A)} 

is a connected, simply connected 2 -dimensional manifold, the same properties hold for the set 

Then it follows that Exp : N 35 —>■ M 35 is a diffeomorphism. □ 
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i 

1 

2 

3 

4 

5 

6 

7 

X 

X 

X 

X 

— X 

— X 

— X 

— X 

y 

-y 

y 

-y 

y 

-y 

y 

-y 


Table 8 Action of e® in the plane {z = 0} 



Fig. 9 Stratification of M' 


4.3 Stratification of the set M' 


Define subsets M' C M', j = 1,..., 40, as follows: 

— For j £ {1, 9,17, 21, 25, 29, 33, 35, 39} , the sets Mj are given by Table 

— For the rest j the sets Mj are given by equalities (4.221-(4.25 I: 


(Mj) = Mj+„ 

i = 1,... ,7, 

e^® 

i= 1,2,3, 

e^+® (Mj) = Mj+„ 

i= 1,2,3, 

(Mj) = Mj+i, 

j = 33,39. 

Lemma 16 A stratification of M' is given as: 

m' 

O II 

□ 

II 


(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 


Proof Follows from Lemma 11 and the description of the action of reflections e® in the plane {« = 0}, see Table 

□ 


m 


Stratification (4.261 is shown in Figure]^ 

Theorem 4.1 For any i = 1,..., 40, the mapping Exp : N[ —> M[ is a dijfeomorphism. 

via the symmetries e® of the exponential mapping. 


Proof Follows from Lemmas 


12 


15 


□ 


Define the following important sets: 
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Fig. 10 Cut Locus 


— the cut locus Cut = {Exp(A, tcut(A)) | A G C} , 

— the first Maxwell set 

Max = {qi G M I 3 minimizers q'{t) ^ t € [0,ti], such that = q"{ti) = qi} . 

— the first conjugate locus Conj = |Exp(A, t™"’’(A)) | A G cl , 

— the rest of the points in M' compared with Cut, i.e., Rest = M'\Cnt. 

We have the following explicit description of these sets: 

Cut = u{M' I i=l,...,34}, 

Max = u{M' I t= 1,...,20,29,...,32}, 

Conj n Cut = u{M' I i = 21,..., 28, 33, 34}, 

Rest = u{M' I t = 35,...,40}, 

Thus we get the following decomposition of the sets M'-. 

M' = Cut U Rest, 

Cut = Max U (Conj n Cut). 

The global structure of the cut locus is shown in Figure From our analysis of the exponential mapping, we 
get the following description of the cut time and the optimal synthesis on SH(2). 

Theorem 4.2 We have the following explicit description of the cut time, tcut(A) = t(A) for any \ £ C. In detail: 

A G Cl tcut(A) = tf""(A) = 4K{k), 

A G Ca tcut(A) = tY^^{X) = 4:kK{k), 

A G C4 Cut (A) = t}°"j(A) = 27r, 

A G C*3 U C 5 => icut(A) = -\-oo. 


-10 
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the case X £ C 4 


Proof If A £ C\C' 4 , then we know from Theorem 3.3 that tcut(A) = t(A) = ti^[X). It remains to consider 


Let A £ C 4 , then qt = Exp(A,t) = (t, 0,0). For any t £ [0, ti], 


the point qt is connected with qo by a unique geodesic Exp(A^,s), s £ (0,si], with (A^ 


t\ = t^) = 27r, 

, si) £ N, namely 


(A^,si) = (A,t) £ A^ 3 g for t £ (0,27r), and (A^,si) = (A,t) £ A ^33 for t = 27r. Thus the geodesic qt, t £ [0,ti] 


IS a mmimizer. 


It follows that tcut(A) = t(A) = t™ 
equality for \£ C\. 


^(A) = 27r for A £ 


Cl 


By applying a reflection e*, we get a similar 

□ 


From the above description of the structure of the exponential mapping, we get the following statement. 

Theorem 4.3 

1. For every point gi £ M U Rest, there exists a unique minimizer q{t), t £ [0,ti], for which the endpoint 
q{ti) = qi is neither a cut point nor a conjugate point. 

2. For any point q\ £ Max, there exist exactly two minimizers that connect qo to q\ for which q\ is a cut point 
but not a conjugate point. 

3. For any point q\ £ Conj n Cut, there exists a unique minimizer that connects qo to q\ for which q\ is both 
a cut and a conjugate point, but not a Maxwell point. 


5 Sub-Riemannian Caustics and Sphere 


In [ 8 ] we presented plots of sub-Riemannian sphere and sub-Riemannian wavefront in the rectifying coordinates 
(Ri, R 2 ,z). Here we perform another graphic study of the essential sub-Riemannian objects, i.e., sub-Riemannian 
caustic and sub-Riemannian sphere. Recall that the sub-Riemannian caustic which is the first conjugate locus is 
given as: 

Conj = |Exp ^A, tj°"'’(A)j I A£C'|. 


The caustic is presented in Figure 11 The component starting at (0,0,0) is the local component of the caustic 


whereas other two parts on right and left side are the parts of the global component of the first caustic. The red 
colored surface inside the local and global components of the caustic is the cut locus whereas we see that the 
boundary of cut locus forms the boundary of the caustic. A zoomed version of the local component of the caustic 


is separately shown in Figure 12 It is evident that it is a four cusp surface as predicted in [?]. A combined plot 


of first and second caustic is also shown in Figure fl3| Note that in the local component of the caustic, the first 
caustic is solid and the second caustic is transparent whereas in the global component of the caustic, the second 
caustic is solid and the first caustic is transparent. The sub-Riemannian sphere S]i{qo',R) at qo is the set of 
end-points of minimizing geodesics of sub-Riemannian length R and starting from qo'. 


Sji = {Exp(A, R)gM 
The following plots are presented: 


A £ (7, tcut(A) > R} = {(J £ M I d{qo,q)=R}. 


1. Sphere of radius R = tt (Figure 141, 

2. Sphere of radius R = 2?! (Figure 15 I, 


3. Intersection of the cut locus with the hemisphere z < 0 of radius R = tt (Figure 16 1 

4. Intersection of the cut locus with the hemisphere z < 0 of radius R = 27r (Figure 

5. Intersection of the cut locus with the hemisphere z < 0 of radius R = Stt (Figure 

6 . Matryoshka of hemispheres z < 0 of radii R = tt and R = 2n (Figure [T^. 


181 


6 Conclusion 

The global optimality analysis and structure of exponential mapping for the sub-Riemannian problem on the 
Lie group SH(2) was considered. We cutout open dense domains by Maxwell strata in the preimage and in 
the image of exponential mapping and prove that restriction of the exponential mapping to these domains is a 
diffeomorphism. This fact leads to the proof that the cut time in the sub-Riemannian problem on the Lie group 
SH(2) is equal to the first Maxwell time. We then describe the global structure of the exponential mapping 
and obtain a stratification of the cut locus in the plane 2 = 0. Consequently, the problem of finding optimal 
trajectories from any initial point qo £ M to another point q\ £ M, 2 yf 0 is reduced to solving a set of algebraic 
equations. Summing up, a complete optimal synthesis for the sub-Riemannian problem on the Lie group SH(2) 
was constructed. 
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Fig. 11 Sub-Riemannian caustic and cut locus 



Fig. 12 Local component of sub-Riemannian caustic and cut locus 



















Cut Locus and Optimal Synthesis in Sub-Riemannian Problem on the Lie Group SH(2) 


29 



Fig. 13 Sub-Riemannian first and second caustic 



Fig. 14 Sub-Riemannian sphere of radius R = tt 









































30 


Yasir Awais Butt, Yuri L. Sachkov, Aamer Iqbal Bhatti 



Fig. 16 Intersection of the cut locus with the hemi¬ 
sphere 2 : < 0 of radius R = tt 


Fig. 17 Intersection of the cut locus with the hemi¬ 
sphere 2 < 0 of radius R = 2tt 




Fig. 19 Matryoshka of hemispheres 2 < 0 of radii R = tt and R = 2n 
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